Superfluidity of bosons on a deformable lattice 
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We study the superfluid properties of a system of interacting bosons on a lattice which, moreover, 
are coupled to the vibrational modes of this lattice, treated here in terms of Einstein phonon model. 
The ground state corresponds to two correlated condensates: that of the bosons and that of the 
phonons. Two competing effects determine the common collective soundwave-like mode with sound 
velocity v, arising from gauge symmetry breaking: i) The sound velocity vq (corresponding to a 
weakly interacting Bose system on a rigid lattice) in the lowest order approximation is reduced due 
to reduction of the repulsive boson-boson interaction, arising from the attractive part of phonon 
mediated interaction in the static limit, ii) the second order correction to the sound velocity is 
enhanced as compared to the one of bosons on a rigid lattice when the the boson-phonon interaction 
is switched on due to the retarded nature of phonon mediated interaction. The overall effect is that 
the sound velocity is practically unaffected by the coupling with phonons, indicating the robustness 
of the superfluid state. The induction of a coherent state in the phonon system, driven by the 
condensation of the bosons could be of experimental significance, permitting spectroscopic detections 
of superfluid properties of the bosons. Our results are based on an extension of the Beliaev - Popov 
formalism for a weakly interacting Bose gas on a rigid lattice to that on a deformable lattice with 
which it interacts. 

PACS numbers: 05.30.Jp, 67.57.Jj, 63.20.Mt, 67.90.-l-z 



I. INTRODUCTION 



The weakly interacting Bose gas has been studied the- 
oretically in great deal over the past fifty years; mainly 
in view of understanding the superfluid properties of "^He 
and its rich phase diagram in the temperature - pressure 
parameter space. These studies provided a qualitative 
description of the low temperature collective soundwave 
like spectrum, the depletion of the condensate, etc. El In 
more recent years attempts have been made to exam- 
ine other systems than "^He which potentially could show 
such superfluid properties. The superfluidity of excitons 
in semiconducting materials such as CU2O, presents one 
of those novel systems.a There, the interaction between 
the excitons with acoustic phonons seems to play a key 
role in establishing such a superfluidity of the bosonic 
excitons. In fact it involves the appearance of a coher- 
ent crystal displacement fipild which enables a moving 
exciton-phonon condensateJj 

In the present work we want to address ourselves to a 
general situation and consider the system in which the 
itinerant bosons, apart from an intrinsic repulsion be- 
tween them, are coupled to some bosonic degrees of free- 
dom, which for brevity we shall term phonons. For in- 
stance, this situation might be realized in short coherence 
length superconductors which are controlled by the fluc- 
tuations of the phase of the order parameter rather than 
of its amplitude. In such systems one expects relatively 
long lived local pairs of electrons which can be considered 
as hard-core Bosons on a lattice and which are coupled 
to the dynamical deformations of this underlying lattice. 



Ion-channelling experiment^ on High- Tc cuprates (show- 
ing a drastic increase in the critical angle for Rutherford 
back-scattering just. below Tc), as well as optical absorp- 
tion measurementsu (showing a substantial increase in 
the phonon intensity just below Tc) indicate that, upon 
entering the superconducting state, the uncorrelated mo- 
tion of the local lattice vibrations might get correlated, 
resulting in the emergence of an acoustic branch in the 
phonon spectrum. 

One also might think about systems, in which spin or 
pseudo-spin (orbital) degrees are coupled to the lattice. 
If, upon lowering the temperature, a long-range order due 
to the continuous symmetry breaking occurs in those sys- 
tems, while the symmetry restoring variable is being cou- 
pled to the phonons, then the system will show the same 
physics associated with macroscopic quantum effects as 
the model we are going to deal with. 

The study presented in the present paper is, to our 
knowledge, the first generalization of the known field- 
theoretical treatment of the weakly interacting dilute 
Bose gas when the effective two-body potential is sup- 
plemented by a time dependent retarded part stemming 
from the phonon mediated retarded interaction. 

We shall study here the feasibility of such a phe- 
nomenon on the basis of a simple model of itinerant 
bosons on a deformable lattice. In the present work we 
consider lattices susceptible of local lattice deformations. 
Such local deformations occur in the systems built up 
of molecular units and manifest themselves as dynam- 
ical deformations which are a priori spatially uncorre- 
lated and are described by the Einstein phonon model. 
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Generalization of such a scenario to globally deformable 
lattices, described by acoustic phonon modes, will be 
briefly mentioned in the concluding part of the paper 
and will be subject of a future publication. In section 2 
we shall define the model and study it within the stan- 
dard Bogoliubov schemes which we extend to two inter- 
acting boson fields: that of the itinerant bosons (the lo- 
cal electron pairs or electron-hole pairs respectively) and 
that of representing the lattice degrees of freedom (the 
optical or acoustic phonons respectively). This approx- 
imation treats the scattering processes among conden- 
sate quasi-particles and between condensate and out-of- 
condensate quasi-particles, but neglects scattering among 
out-of-condensate quasi-particles. These latter processes 
are important if the system are not very dilute and if 
the temperature is finite. A consistent theory which can 
satisfactorily deal with, that situation is the so called 
Beliaev-Popov theoryQ'Q, which, in section 3, we shall 
generalize from a weakly interacting Bose gas on a rigid 
lattice to one where also the interaction between the 
bosons and the lattice vibrations are taken into account. 
In section 4 we address the question of the renormal- 
ization of the sound velocity of the collective modes in 
the hydrodynamic regime at temperature zero. We shall 
illustrate the subtle compensation between two compet- 
ing mechanisms arising from the boson-phonon coupling 
which while leading to an enhancement of the boson mass 
it also leads to an increased rigidity of the phase of the 
condensate. As a result the sound velocity is practically 
the same as that for a weakly interacting Bose gas on 
a rigid lattice. Concluding remarks and an outlook on 
further studies are presented in section 5. 



II. THE MODEL AND ITS SOLUTION WITHIN 
A BOGOLIUBOV SCHEME 

We consider a system of bosons on a lattice having a 
tight binding spectrum £q = zt[l ~ l/'i{cosqx + cos qy + 
cos^z)] (From now on we shall set zt — 1, measuring 
all energies in units of the half bandwidth) which in the 
long wavelength limit is given by Eq — /2M with M 
being the the bare boson mass. The interaction between 
bosons is characterized by a coupling constant g and that 
between the bosons and the phonons by awo, i^o denoting 
the frequency of the optical phonon mode. The Hamil- 
tonian for such a system is then given by 

H = Y^ eq5jj6q -I- ujo '^{a^a^ + i) 
q q 

k,k',p 
aujQ , t , r t 



(1) 



k,q 



where b^"^ and Oq ^ denote the boson and phonon annihi- 
lation (respectively creation) operators, Cq = £q — /i and 



fj, being the chemical potential. Let us now assume the 
existence of a condensed state not only for the bosons but 
also for the phonons and consequently make the Ansatz: 



6k = bu 



bS, 



k.O, ak — Ok 



0(5, 



k.O 



(2) 



By substituting expressions (g) into the Hamiltonian (|i|) 
and requiring that the terms linear in 6k and Oq vanish 
we find the following relations 

ab'^ r— 
a = — = — a\ Nrin 



N 

N 

fj. guc, g = 9 
finally yielding the following Hamiltonian 



2a^wo 



H = Hn 



H 



+ H 



B-B 



(3) 



(4) 



and where 



-^0 = X! <^q^q^q 



Wo 



6q6-q 



q 



with 



2gnc — 2a ajQUc = £q + gric. 



(5) 



(6) 



The Hamiltonian Hq describes the Bogoliubov quasi- 
particles which are hybridized with the phonons. In 
the condensed state the single particle excitations hy- 
bridize with the density fluctuations due to the broken 
gauge symmetry and give rise to the collective excita- 
tion spectrum of the bosons with a sound-wave like spec- 
trum. The phonons which initially couple to those den- 
sity fluctuations hence get hybridized with those collec- 
tive modes. The second term in the Hamiltonian Eq. (^) 
describes that coupling between the density of the out- 
of-condensate particles with the phonons. 



Hb-p — - 



k^k+q [a^ + a_q] 



blh 



(7) 



The last term in the Hamiltonian Eq. (j^, Hb-b, de- 
scribes the interaction between the out-of-condensate 
particles. 

Let us next define a generalization of the standard Bo- 
goliubov transformation 

aq = Ul,q6q + U2,qaq + Wl,q6t.q -|- W2,qa-q 
/3q = Ul^J>q + ?22,qaq + Wl^qSlq -|- W2.qa^q 
alq = Wl,q6Lq + M2,qa-q + Wl,q6q + W2.qaq 

W2,qa^q + Wl.qOq + W2,qaq. (8) 
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With the requirement that the new operators satisfy the 
standard canonical commutation relations for Bosons, we 
obtain the following relations between the coefRcients of 
this transformation 



i,q 



^2,q 



"l,q 



2 



"l,q "'2,q ~ "l.q ~ '"2,q 

"l.qMl.q + M2,qM2,q " I'l.qWl.q " ■f^2,q'l'2,q 
'"l,q'l'l,q + U2,qW2,q " Wl,qWl,q " 1'2,qW2,q 

The transformation, Eq. (^), diagonalizes the Hamilto- 
nian _ffo, bringing it into the form 



(9) 



Hq^^ (t^l.qftqaq + W2,q/3^/?q) 



(10) 



The eigen-frequencies wi.q and Ci;2.q as well as the co- 
herence factors are determined in the usual way by 
diagonalizing the secular equation for the variables 
Qfq, aLq, /3q, /3ln which Icads to: 



UJ 



l,2,q 



^Iq 



_ CJl,q ± £q 
S-^t^l.q^q 

_ ^l,q ± ^0 
2-^(^1, qti^O 

'^2,q ± 



,212 



IGa^WgCqric 



"l:q 



-'2,q- 



'l,q 



"l:q 



"2,q 



'l,q 



0,q 



"2,q 



2^^:;^;^ V '^2,q ~ '^ii 



W2,q±Wo /'^2,q- 



^Iq 



2^(x'2,qt^0 y t^2,q ^ ' 



(11) 

(12) 
(13) 
(14) 

(15) 



^l,q 



and where — \/ £q[£q + ^gn^ denotes the spectrum 
of Bogoliubov quasi-particlcs in the absence of coupling 
with phonons. 




(7t,Jt,Jl) 



FIG. 1. Sketch of the excitation spectra ii'i,2,q, Eq. 
in the lowest order approximation (solid line). The excita- 
tion spectra in the absence of boson-phonon coupling is rep- 
resented by dashed lines. 

In Fig. 1 we illustrate the two branches of eigen- 
frequencies wi(2),q. When -Esq is not close to the phonon 
frequency wq, then the two normal modes consist of a Bo- 
goliubov excitation and a phonon mode; one of the modes 
being mainly a phonon and the other a Bogoliubov ex- 
citation. For i?Bq close to wq, neither mode is predomi- 
nantly a phonon or a bosonic excitation. The mode which 
is a phonon (Bogoliubov excitation) for Ebo^ ^ be- 
comes a Bogoliubov excitation (phonon) for ^ ^-Sq- 
In the long-wavelength limit the lower branch reduces to 
a sound-like dispersion 



t^i.q - wog, q 







(16) 



with a characteristic sound velocity 



(17) 



which is reduced in comparison to that where the boson- 
phonon coupling is absent. This reduction is determined 
by the static part of the attractive phonon mediated in- 
teraction. 

The set of relations, Eqs. (11 - 15) represent the lowest 
order approximation to our problem, i.e., it takes into ac- 
count collisions between condensate quasi-particles, be- 
tween condensate and out-of-condensate quasi-particles, 
but totally neglects scattering among out-of-condensate 
quasi-particles as well as scattering out-of-condensate 
quasi-particles on phonons. As the concentration of 
quasi-particles increases these latter become important. 
We therefore shall have to include in our study the ef- 
fect of Hb~p and Hb-b- In the next section, wa.ahall 
generalize the Beliaev- Popov approximation (BPA)lIi3 for 
the case when the bosons, coupled together by the two- 
body potential, are moreover coupled by the phonon- 
mediated retarded interaction. Since in the BPA the 
second-order self-energies are built out of the mean-field 
propagators we shall define all the Green's functions in- 
volved in our perturbation theory and give their expres- 
sion within the mean-field approximation. These propa- 
gators can be defined as the elements of a matrix Green's 
function composed of a four-component bosonic field 



q' 



7 h'^ 

^qj q; 



^}, the components of which are how- 
ever related among each-other by symmetry relations. 
We define the canonical 2x2 matrix Green's functions for 
bosons, with diagonal (Q) and non-diagonal (Q) compo- 
nents; the phonon Green's functions (I?) and the phonon- 
boson Green's function (7i), that describes the hybridiza- 
tion of these two bosonic excitations in the condensed 
state 



3 



2?q(w) = ((flq + al^\al + a_q))^ , (18) 

^-Lc^iuJ) = {{bc^\al + a-c^))u ■ 

The mean-field Green's functions, introduced above, are 
obtained from the diagonalized Hamiltonian, Eq. (10), 
with the help of the inverse canonical transformation 



Wl,q/3q 



flq = U2,qQ;q + ?l2,q/3q - ^^2,qa!_q " W2,q/3lq 
^Lq = -I'l.qftq " Wl,q/5q + Wl.qttLq + "l.q/^lq 
aLq = -1'2,qaq " W2.q/3q + W2,qaLq + ?l2,q/3lq. 



(19) 



This leads to the following expressions for those mean- 
field Green's functions 



i,q 



t'q.t:; = ""l.q^'l^q 2 _ , ,2 1" i^l.q' " 



*,t^2,q,U,'!;J , 



Wq.i^ = ('«2,q - 'y2,q) 



2c^l,q 



■"l,q 



Wl^q 



_W - Wl^q 
-I- [a;i,q,U,?; ^,tJ2,q,M, W] 



Ll! + Ul , 



(20) 



In Fig. 2 we express graphically the phonon "D and 
the boson-phonon TC Green's functions in terms of G, the 
normal and G, the anomalous component of the boson 
Green's functions. The choice of this graphical repre- 
sentation is convenient when constructing all the second- 
order self-energy diagrams, as will be discussed in the 
next section. The full straight lines with two arrows 
in the same (opposite) direction represent the normal 
(anomalous) component of boson Green's function. The 
thin wavy and bold wavy lines stand for the bare and 
the mean field Green's functions of the phonons. The 
boson-phonon Green's functions are built up of a wavy 
line and an outgoing straight line. The dot represents the 
boson-phonon vertex and the zig-zag line the condensate 

1 /2 

(given by a factor ric )■ 



/\/\/\/\ = /\y\/\/\ + 



_^\/\/\/\ + 



FIG. 2. Graphical representation for the mean field Green's 
functions "D and Ti,. 



III. GENERALIZED BELIAEV-POPOV THEORY 

A consistent theory which goes beyond the mean-field 
approximation and takes into account collisions between 



out-of-condensate quasi-particles waSpdeveloped by Be- 
liaev for a dilute Bose gas at T = O.Q A generalisation 
of this for finite temperature was done by Popov .0 The 
Beliaev-Popov approximation is the next step beyond the 
lowest order gapless theory and provides a consistent 
second-order self-energy approximation for the weakly 
interacting dilute Bose gas. It gives a gapless excita- 
tion spectrum with the velocity equal to the macroscopic 
speed of the sound at T = 0, in agreement with the find- 
ings by Gavoret and NozieresO which is valid to all orders 
in a perturbation theory treatment. For details on the 
underling physics of the BEA we refer the reader to the 
review by Shi and Griffin. E2l 



A. Effective interaction and higher-order 
contributions the self-energy 



In this section the BPA is generalized to the case where 
the quasi-particles, interacting via a standard two-body 
potential, are moreover coupled by a phonon mediated 
retarded interaction. To construct the second-order self- 
energy diagrams we define an effective interaction as the 
sum of a two-body t-matrix and a phonon mediated re- 
tarded interaction a'^ujQ'Do{uj), with Vq^lu) being given 
by the bare phonon Green's function. 



FIG. 3. Graphical representation of the effective bo- 
son-boson interaction. 

The effective interaction is illustrated in Fig. ^ as a 
hatched line. The dashed line stands for the bare in- 
teraction potential and the non-hatched interaction line 
denotes the two-body i-matrix. 

In Fig. 4 we reproduce all the Beliaev-Popov diagrams 
for the noHHal (S) and anomalous (E) components of 
self-energy.ll3 Each interaction line (hatched line) can be 
either a t-matrix or a phonon propagator. The diagrams 
are grouped in a way which is different from the con- 
ventional one. As it turns out this type of grouping is 
more convenient when it comes to summing up in a com- 
pact form the subset of self-energy diagrams involving the 
phonon mediated interaction. We use the standard nota- 
tion of Green's function lines having cuts which indicate 
diagrams which are already included in the t-matrix and 
should be not counted again. The mean field approx- 
imation, discussed in the previous section, retains only 
the first two diagrams from the group al and the first 
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diagram from the group a2 in the normal and anoma- 
lous components of the self-energy respectively. The t- 
matrix is then being approximated by the s-wave scat- 
tering length. 

Each second-order self-energy diagram consists of three 
contributions. The first contribution arises when both 
effective interaction lines are replaced by the i-matrix. 



They have been extensively discussed in Ref. The 
only difference between that and the case we want to 
study here is the four pole structure of the mean field 
Green's functions [see Eq. (^)]. The second contribu- 
tion involves both the <-matrix and the phonon mediated 
interaction. Finally, the third contribution is exclusively 
due to the phonon mediated interaction. 



al) 




i 



a2) 







^^'I'^^^'j^ ^^'l^^'j^ ""l^"^ "^i^"^"!^ ^^^j^--<-|^ 



tX ±± ±± XX xl il 



ci) ^/^fT^ ^ ^jj^ 



dl) 



^^^^ 



c2) 



d2) 



^^^^ 
^ > < ^^^^^^^ ^^^ 



► Ig) < ^^^^^^ ^^^^ 



FIG. 4. Graphical representation, up to second-order in the effective interaction of all contributions to the normal and 
anomalous self-energy. The hatched lines represents the effective interaction shown in Fig. 3. 



To sum up in the more compact form the subset of 
diagrams involving the boson-phonon coupling we use 
the graphical representation for the mean-field Green's 
functions, derived in the previous section [see Fig. 2]. 
We focus our discussion on the third subset of diagrams 
defined above. The second subset of diagrams, involving 
both types of interactions, are summed up using the same 
procedure, ft turns out that the summation of the first 



four diagrams from group bl together with the fourth 
diagram from group al [see Fig. 4] is equivalent to re- 
placing the bare phonon propagator in the last diagram 
by the mean-field propagator. As a result one arrives at 
the first diagram shown in Fig. 5. Similarly, the sum- 
mation of the first four diagrams from the group cf in 
Fig. 4 leads to the second diagram in Fig. 5. The latter 
one represents the convolution of the Green's functions 
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describing the boson-phonon hybridization. The other 
diagrams in Fig. 5 and the ones for the anomalous self- 
energy, represented in Fig. 6, are obtained in the similar 
way. 




FIG. 5. Graphical representation, up to second-order in 
the phonon mediated boson-boson interaction, of the normal 
self-energy E due to the phonon mediated interaction. 



4% 




FIG. 6. Graphical representation up to second-order in the 
phonon mediated boson-boson interaction, of the anomalous 
self-energy E due to the phonon-mediated interaction. 

We would like to point out that, in the present treat- 
ment, the renormalization of the phonon mode as well 
as of all relevant vertex corrections is consistently taken 
into account up to the second-order in the effective in- 
teraction. The phonon mode is indeed renormalized by 
the density fluctuation of bosons. In the condensed state 
the density response function consists of two contribu- 
tions: a quasi-particle part and a part describing density 
fluctuation of out-of-condensate particles. The phonon 
renormalization due to the first contribution is already 
present in the extended Bogoliubov scheme, leading to 
the renormalized (mean-field ) propagator of phonons. 
By replacing the effective interaction by phonon propaga- 
tors, one verifies that the first diagram from group al in 
Fig. 4 together with last two diagrams from group bl de- 
scribe the phonon mode renormalization arising from the 
density fluctuations of out-of-condensate particles. They 
are described given by the "particle-hole" bubble dia- 
grams built up from normal and anomalous mean-field 
(Bogoliubov) Green's functions of the bosons. One far- 
ther verifies, that the diagrams from group dl together 
with diagrams of identical topology from group cl from 
Fig. 4 generates the vertex corrections to the boson- 
phonon as well as to the boson-boson vertex. For ex- 
ample, when all the effective interaction lines in these 
diagrams are replaced by the phonon mediated interac- 
tion, one obtains the vertex corrections to the phonon- 
boson interaction. Similarly appear the corrections to 
the boson-phonon vertex due to the boson-boson inter- 



action (these diagrams are generated by replacing one 
effective interaction line by the boson-boson interaction 
line). When all the phonon lines in the above mentioned 
class of diagrams are replaced by boson-boson interaction 
lines one obtains the renormalization of the boson-boson 
interaction. In a similarly way we can generate the dia- 
grams renormalizing the boson-boson interaction by the 
phonon-boson interaction. 

The normal and anomalous components of the Bo- 
son Green's functions are expressed in terms of the self- 
energy via the Dyson-Beliaev equation 



Gq.i 



(21) 



where Cq is the bare excitation spectrum measured from 
the chemical potential ^. In Eq. ( |l|) the following no- 
tations are introduced: 



\lu-A. 



(22) 



Here Aq^u) is an antisymmetric function in lu, while Sq^i^ 
and Sq.i^ are symmetric functions of uj. Invoking the 
Hugenholtz-Pines theorem (see the following sub-section) 
we have the relation /j, = Eo,o — £o,o — •S'o.o- 

The existence of the two branches in the excitation 
spectrum, Eq. ( pd]) leads to three different contribu- 
tions to the self-energy, describing the intra-band and 
inter-band scattering processes. Based on the mean-field 
Green's functions, Eq. (pO|), we evaluate the above in- 
troduced dynamical quantities at zero temperature and 
arrive at the following result: 



A. 



^a,/3(k,q)t^ 



,k + t^/3,k+q)^ 



5 =1 V 



ScfsO^, q)('JJa,k + t^/3,k+q) 
^2 - (LUa,k + t^/3,k+q)^ 

J_ \p Mg^pjk, q){uJaM + ^/3,k-|-q) 



k,a,/3 



(23) 



a {(3) denotes the branch index, f3 < a — 1,2, the cor- 
responding vertices Aa^f3, Sa,i3, and Ma^p which are ex- 
pressed in terms of the coupling constants and coherence 
factors and are given in Appendix A. 

The renormalized quasi-particle spectrum is given by 
the poles of the dressed boson Green's function, or equiv- 
alently by the zeros of Dq^^^ . In the next section we dis- 
cuss the long- wave length excitation spectrum of the sys- 
tem, focusing on the effect of boson-phonon coupling on 
the sound- velocity renormalization. 
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B. The Hugenholtz-Pines theorem 

One of the fundamental results in thej-theory of Bose 
systems is the Hugenholtz-Pines theorenJll which relates 
the chemical potential to the self-energy in the condensed 
state via the relation 



(24) 



This theorem implies that the excitation spectrum of the 
system is gapless. However in certain approximations, 
for example the Hartree-Fock-Bogoliubov approximation 
(HFB), higher-order terms in the chemical potential are 
retained, while being neglecting in the self-energies. Such 
approximations are not "consistent" to any given order 
in the coupling constants and can generate an unphysical 
gap in the excitation spectrum of the system. Let us now 
verify this point for our treatment. We first calculate the 
chemical potential from the Heisenberg equations of mo- 
tion for the Bose field (analogous to the Gross-Pitacvskii 
equation) and compare the result to that obtained from 
the Hugenholtz-Pines theorem. Our further treatment 
closely follows that of Ref. 
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We start from the exact Heisenberg equation of motion 
for the Bose field b„ 



ibr, 



^k+p^k+qfop 



k.p 



X]^k+q[ak + a-k]. 



(25) 



k,q 



The next step is to separate out the condensate part 
as in Eq. (|) and treat the boson interaction term 
in a self-consistent mean-field approximation (HFB 
approximation) .Ej As a result we arrive at the following 
equation 



^5 = .gn,6+^^[2(6^6q) + (StS, 



q 



qj/ J 



(26) 



where Ir /N = ric is the condensate fraction and n is the 
boson density. The above equation yields the following 
result for the chemical potential: 

fi — fii + ^2, Ml — 9^1 = g[n + m] — aujoh , (27) 

ri, ifi denote the density and anomalous density of out-of- 
condensate particles, respectively and h is the strength of 
boson-phonon hybridization. These quantities are given 

by 



^q/ I 



^ = ;7=]^E(^[«q + «-q]) 



(28) 



So far we did not make any approximation for the 
boson-phonon interaction. We next show that if all the 
above defined averages are calculated self-consistently 
from the corresponding mean field Green's functions, 
then the derived chemical potential Eq. (27) will coincide 
with that obtained from the Hugenholtz-Pines theorem. 

In Fig. 7 we represent graphically the contributions 
to the chemical potential, Eq. ([27|), arising from the cou- 
pling of the bosons with the phonons. Putting this con- 
tribution together with the first contribution to /i2 in 
Eq. (27) we obtain /i2 = A -I- A. A and A are the con- 
tributions corresponding to the fourth diagram from the 
group al and the second diagram from the group o2, 
presented in Fig. 4. 



n 



FIG. 7. Graphical representation of second-order contri- 
butions to the chemical potential due to the coupling with 
phonons. 

Now we calculate the chemical potential from the 
Hugenholtz-Pines theorem. Examining the various dia- 
grams of I]o,o and Sq^o we find that most of them cancel 
pairwises, giving the following result 



A*HP — gn + g[n + m^] + A — A + B 
gric \ - _l_ 
^ 977 ■ 





Tor — TO 



N 



(29) 



Tor denotes the " renormalized" anomalous densityS and 

B - ^ E [<^Mg - «'^o2'o(^)]'{a^a-^ - \QJ\^} . 

q J 

Using the Dyson-Beliaev equation ( pT| ) we verify the fol- 
lowing identities 

^ = ^ E / '^^[S - o?^lT^^{^)\Qu. = 2i . (30) 

q 

Substituting the expressions, Eq. ( |30| ) into Eq. ( |29| ) 
we see that the chemical potential derived from the 
Hugenholtz-Pines theorem coincides with the one ob- 
tained from the Heisenberg equations of motion, Eq. ([2^). 

IV. RESULTS AND DISCUSSION 



In the preceding section we have given the expressions 
for the self-energies in the Beliaev-Popov approximation. 
We now use these results to calculate the second-order 
corrections to the sound velocity and the chemical po- 
tential. 
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A. Depletion of the ground state 

We first discuss, the depletion of the ground state. 
From the Green's function Eq. ( ]20| ) we determine the 
density of out-of-condensate particles n, and thus obtain 
the self-consistent equation for the density of particles in 
the condensate: 



k 



l.k 



(31) 



f [gn-c 



c] [wo + £^k 



2^k(wi,k + W2,k) 



+ 



Wi,k + '^2,k 



Eq — \/sq[eq + 2gnc] is the spectrum of Bogoliubov 
quasi-particles with reduced scattering length g [see 
Eq. (|^)]. The self-consistent equation for the conden- 
sate fraction, Eq. is solved numerically. In Fig. 8 
the relative depletion r — n/n as a function of density 
is presented for various values of the boson-phonon cou- 
pling constant a. As clearly evident from this figure, for 
vanishing a, the relative depletion follows the canonical 
behavior, which is characteristic for the dilute weakly 
interacting Bose gas. For a non-zero coupling constant 
however, r remains finite in the limit n — )■ 0, indicating 
a linear in density contribution to the depletion of the 
ground state. 
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FIG. 8. Relative depletion of the ground state r — n/n as 
a function of density for different values of coupling constants 
{a =0, 1, 1.5, and 2 - solid line, dotted, dashed and long 
dashed lines, respectively) for ujq — 0.2 and g=3. 

For 71 ^ 1, expanding the denominator in the last sum 
of Eq. (31) and keeping the leading order contributions 
in Uc we arrive at the following expression 



1 



N ^ 

k 



^k — ^^k 



(32) 



(cjo + £^k)2. 

The first term in Eq. (Isl) gives rise to the canonical ri^J"^ 



contribution to the ground state depletion while the sec- 
ond term leads to the linear in density contribution, and 
is due to the boson-phonon coupling. 



B. The chemical potential 

The second order contribution to the chemical poten- 
tial is given as a sum of two contributions: /J.2 = l^b + f^p- 
The first one is due to the boson-boson interaction and 
the second one stems from the boson-phonon interaction: 



fib = gin + ?%,) 
fip = —auoh = —a~uj, 



—T 

2N ^ 

k 

E 



£k ^0 



2, ,2 



E\^ UJl,k + 
gk 

-Ek(Wi,k + t^2,k) ' 



gn-c 



(33) 



In the normal state, where Uc = 0, the second contribu- 
tion fip coincides with the density independent negative 
shift of the chemical potential due to the coupling with 
phonons. However, in the condensed state this contribu- 
tion becomes density dependent, via the density depen- 
dence of the spectrum in the condensed state. The den- 
sity dependent part of this contribution increases with in- 
creasing density and hence gives rise to an increase of the 
system's compressibility. In Fig. 9 we plot the density 
dependent part of this contribution (fip — fJLp{n) — /io(0)) 
for several coupling constants. 
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a. 




0.04 



0.02 



FIG. 9. Second order contribution to the chemical poten- 
tial due to the boson-phonon coupling (/ip) versus density (n). 
[a = 1, 1.5, and 2 - dotted, dashed and long dashed lines, 
respectively; too = 0.2 and g=3. 

The total second order correction to the chemical po- 
tential is presented in Fig. 10 which also shows its in- 
crease with a. We would like to emphasize that for a = 
and for n ^ the second order correction to the chemi- 
cal potential shows the canonical n^/^ behavior, while for 
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large a the slope of /i2 is almost linear [see inset on Fig. 
10]. 




FIG. 10. Second order contribution to the chemical poten- 
tial ^2 vs density, (a = 0, 1, 1.5, and 2 - solid, dotted, dashed 
and long dashed lines, respectively, u)^ = 0.2 and g = 3. 




FIG. 11. Chemical potential [i vs density n. (a =0, 1, 1.5, 
and 2 - solid, dotted, dashed and long dashed lines, respec- 
tively, LOQ = 0.2 and g = 3). 

Keeping only terms in leading order in the density, the 
expression for Eq. becomes 



9_ 

N 



E 



^0 + 2^k 







a- 





(34) 



It is the contribution in square brackets in Eq. (Q), which 
leads to the canonical n^^^ behavior of the chemical po- 
tential. On the contrary, the second term, which is due 



to boson-phonon coupling, is linear in density. Upon in- 
creasing a this latter contribution to the chemical po- 
tential increases while the term proportional to n^^^ de- 
creases, because of the reduction of the effective scatter- 
ing length g. This behavior is clearly seen from Fig. 11, 
where the chemical potential = /ii + is illustrated as 
a function of density. For small n the contribution linear 
in density is dominating and the reduction of the chem- 
ical potential due to the phonon-mediated effective at- 
traction in lowest order is almost canceled by the second 
order contribution. As a result, the slope of /x is almost 
independent of a at small densities [see the inset in Fig. 
11], while at higher values of n the chemical potential 
decreases with increasing of a, because of a decrease of 
the contribution proportional to n^^^. 



C. The renormalized sound velocity 

The renormalized quasi-particle spectrum, arising from 
the second order self-energy corrections, are given by the 
poles of the Boson Green's function (|2^), or equivalently 
by the zeros of Dq ( p3| ) . Expanding the self-energies in 
the long-wavelength and low frequency limit around the 
mean field spectrum, we obtain the renormalized sound 
velocity v as 



v = vo[l + X]. 



(35) 



vq = \/ gric/M denotes the mean field sound velocity and 
A the velocity renormalization factor given by 



A = 5/ {2gnc) + gUcS - a. 



(36) 



V(2) 
^0,0 



denotes the second order correction to the anoma- 
lous self-energy. The symmetric and antisymmetric self- 
energies, to leading order in to and q [see Eq. (23)], are 
given by 



-'l,q 



(37) 



Each of the three contributions entering in sound ve- 
locity renormalization factor, Eq. (|36|), contains infrared 
divergences, which are inherent to the Bose condensed 
state. However, as in the case of the standard weakly 
interacting Bose gas, these divergent terms are exactly 
canceled out in the final result for the low-frequency ex- 
citation spectrum. We have checked this cancellation an- 
alytically, picking up the divergent terms both due to the 
coherence factors and the existence of the gapless mode. 

In Fig. 12 the sound velocity renormalization factor 
is presented for various values of the coupling constant 
a. In the absence of boson-phonon coupling A follows 
the canonical square root behavior in the density.t2l An 
increase of a results in an overall enhancement of A and 
moreover leads to a finite value of A, even for vanish- 
ing density. As a result, the renormalized sound velocity 
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[see Fig. 13] is almost unaffected by tfie coupling of the 
bosons with phonons, in contrast to its mean field result, 
which is also presented in Fig. 13. 



<< 1 




FIG. 12. Velocity renormalization factor A versus density 
n. {a —Q, 1, 1.5, and 2 - solid, dotted, dashed and long 
dashed lines, respectively, uio = 0.2 and <; = 3). 




FIG. 13. Renormalized sound velocity v (upper curves), 
mean-field sound velocity (lower curves) versus the density n 
{a =0, 1, 1.5, and 2 - solid, dotted, dashed and long dashed 
lines, respectively, uio ~ 0.2 and g — 3). 

The interesting question now is, why the boson dress- 
ing effects, leading to its mass enhancement, does not 
show up in the velocity renormalization factor. Respon- 
sible for the mass enhancement is the first diagram in 
Fig. 5 for the normal self-energy. It is frequency depen- 
dent and is finite in the limit n — > 0. It gives a positive 
contribution to a and hence a negative contribution to 
the sound velocity [see Eq. (|3^)]. It turns out that in 
the condensed state there are two effects that overcom- 



pensate this negative contribution. First, the effective 
boson-phonon coupling is no longer local and becomes 
momentum dependent because of the coherence factors. 
It is decreased in the condensed state. For the momen- 
tum transfer equal to the wavevector at which the level 
crossing of the bare modes occurs, this effective interac- 
tion is exactly zero. The second main important contri- 
bution competing with this diagram is contained in the 
first diagram illustrated in Fig. 6 for the anomalous self- 
energy. This diagram is positive, thus leading to a posi- 
tive contribution to the anomalous part of inter-particle 
potential, and hence increases the sound velocity. 

In Fig. 14 we present these two contribution for var- 
ious a, rescaling them by a factor {o^ujq)^^. As can be 
seen from this figure there is almost an exact cancellation 
of these two contributions for different a. That explains 
the absence of the density independent negative contri- 
bution to the sound velocity originating from the boson 
mass enhancement. 




-0.02 



-0.04 



FIG. 14. The negative (positive) contribution to the sound 
velocity renormalization factor from normal (anomalous) sec- 
ond order diagrams illustrated by the first diagram in Fig. 5 
(Fig. 6) {a — 1, 1.5, and 2 - dotted, dashed and long dashed 
lines, respectively, cjq ~ 0.2 and <; = 3). 



V. CONCLUSION 

In the present paper we have considered the problem of 
a system of bosons on a deformable lattice which, apart 
from an intrinsic repulsion between them, are coupled to 
local lattice deformations - treated in terms of Einstein 
phonon modes. This leads to a retarded time dependent 
contribution to the effective two-body vertex. When the 
bosons condense, the phonons, being initially coupled to 
the boson density (symmetry restoring variable in the 
case of gauge symmetry breaking), get hybridized with 
the Goldstone mode. The generalized Bogoliubov trans- 
formation has been used to diagonalize the bilinear part 
of Hamiltonian describing the excitation spectra of the 
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system in the lowest order approximation. As a result 
two dispersion branches have been obtained: one describ- 
ing a gapless sound-wave like mode and one exhibiting a 
gap. The two normal modes describe BogoHubov type 
excitations and Einstein phonons. For momenta close to 
where the level crossing of bare excitation spectra occurs 
neither mode is predominantly a phonon or a Bogoliubov 
quasi-particle. 

The ground state depletion is shown to consist of two 
contributions: the first one showing the canonical v?^"^ 
behavior in the density and the second one, which ap- 
pears because of the retarded phonon mediated interac- 
tion, being linear in density. The relative depletion of the 
ground state remains finite in the limit 71 — * because of 
retardation effects. 

The sound velocity obtained in this approximation is 
reduced due to reduction of the repulsive boson-boson 
interaction, arising from the attractive part of phonon 
mediated interaction in the static limit. 

Considering the second order corrections to the chemi- 
cal potential and excitation spectrum, the Beliaev-Popov 
theory has been generalized to our case. For that pur- 
pose an effective interaction has been introduced and 
all second-order self-energy diagrams were constructed 
in terms of this interaction. The short range bare boson- 
boson two-body potential has been renormalized in the 
standard way by summing up the ladder diagrams and 
introducing the two-body ^-matrix with the characteris- 
tic s-wave scattering length. Since the phonon mediated 
boson-boson interaction, vanishes at high frequencies as 
l/w^, this interaction does not require any special renor- 
malization in order to avoid ultraviolet divergences in the 
case of continuum model. Due to the pole structure of 
phonon mediated interaction, the density of bosons no 
longer enters in the expansion parameter and one has to 
assume a small value of boson-phonon coupling constant 
in order to treat the problem by perturbation theory. 

Unlike the lowest order results, the second order con- 
tribution to both the chemical potential and the sound 
velocity shows an increase with the boson-phonon cou- 
pling constant and contains a term linear in the density 
which is exclusively due to the retarded nature of phonon 
mediated interaction. As a result the total chemical po- 
tential and the sound velocity obtained within this theory 
are practically unaffected by the coupling to the phonons. 
This effect is more pronounced in the low density limit 
n ^ 1. One might ascribe this effect to the robustness 
of the superfluid state. 

In the present paper we have restricted oursclf to the 
case when the boson system is coupled to optical phonon 



modes of the lattice. The formalism developed here can 
be easily extended to the case of acoustic phonons. In 
such a scenario the two normal modes correspond to 
two types of sound modes: the acoustic phonon mode 
of the lattice and the Goldstone mode of condensed Bose 
system. The effect of the coupling between these two 
sound modes as well as as well as damping effects can be 
treated within the present formalism. However, in that 
case the coupling with g = phonon mode (representing 
translation of the whole crystal) should be discarded. It 
amounts to neglecting the Hartree type contribution to 
the boson self-energy due to the phonon mediated inter- 
action. 

The present study induces us to speculate on the possi- 
bility of an insulator to superconducting transition when 
the boson-phonon coupling constant is strong. In that 
case the mass renormalization of the bosons varies expo- 
nentially with the coupling constant as exp(Q;^) when we 
consider the normal state of the system. In the superfluid 
state we expect again a phase stiffness of the condensate 
practically unaffected by the coupling of the bosons to 
the lattice vibrations. Possibly such features exist in the 
cuprate superconductors which show a resistivity which, 
upon lowering the temperature, tends to an insulating 
behavior before abruptly dropping to zero when the sys- 
tem becomes superconducting. Similarly, upon entering 
the superconducting state, out of the broad incoherent 
contribution to the angle resolved photo-emission spec- 
tra evolves a sharp resonance peak which could suggest 
well defined quasi-particles in the superconducting state 
and totally diffusive modes in the normal state. Within 
the physics developed in the present study we can spec- 
ulate on an undressing of the bosons as the temperature 
is lowered and the transition from the normal into the 
superfluid state takes place. These and related questions 
will be addressed in some future work. 
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APPENDIX: EFFECTIVE VERTICES 

As we have already mentioned in the main text, each self-energy contribution introduced in Eq. (|2^) consist of 
three different parts describing the intra-band (denoted as 11(22) for the scattering within the lower (upper) branch) 
and inter-band (denoted as 12) scattering processes. The vertices for inter-band scattering are given by 
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^12 (k, q) = Ag'^ric (w?,kWl,k+q - ^^l,k^^l,k+q) - 2wi,kWl,k(M?,k+q - ^^l,k+q) - 2ui,k+qUl,k+q(w?_k " ^^l,k) 
- Agau>o^/n^ (U2,k - V2,k){(wi,k - ■f^l,k)(M?,k+q - '^l.k+q) - ('"l,k + 1'l,k)wi,k+qWl,k+q} 

+ (W2,k+q - W2,k+q){(Sl,k+q - Wl,k+q)(Ml,k ^ '^'l.k) ^ (^l,k+q + Wl,k+q)wi,k?^l,k} 



'LOOHc 



Ul,k'i^l,k(u?,k+q - ^^l,k+q) + ^^l,k+qWl,k+q(w?,k " «?,k) 
2a^Ull^/n^ {(W2,k - V2,k)('"l,k+q + ^^l,k+q) + (■W2,k+q - ^^2,k+q)('«l,k + Vl,k)} {wi,k^^l,k+q + "l.k+qVl.k} 
(U2,k - t^2,k)^(uf^k+q - ^l,k+q) + (W2,k+q - ^^2,k+q)^(w?_k " ^^l,k) 



, 2 2 



+ 2(W2,k - f2,k)(W2,k+q - f2,k+q)(wi,kWl,k+q - Vl,kWl,k+q) 

r 1 ^ 

S'i2(k,q) = 45^nc ui,kUi,k+q - ?^i,k^^i,k+q 

- 4gau>o^/n^ (wi.kUl.k+q - l^l,kt^l,k+q){(W2,k - V2,k)('Ml,k+q + ^^l,k+q) + (■W2,k+q - t^2,k+q)('Ml,k + 1^1,k)}j 
+ a^cul (U2,k - V2,k)^(wi,k+q + ^^l,k+q)^ + ('W2,k+q - ^^2,k+q)^('"l,k + t^l,k)^ 

Ml2(k, q) = 85 V (W?,k^'l,k+q + ^*l,k+q^l,k) - ^^l,k^^l,k(w?,k+q + ^^l,k+q) 

- 2Ml,k+ql'l,k+q(Ml,k + ^'l.k) + 3Ml,k^'l,kUl,k+q^'l,k+q 

+ ^gOLUQ^c 2{(Ui,k - fl,k)(W2,k - t'2,k)wi,k+qt^l,k+q + (wi,k+q - t^l,k+q) (W2,k+q - W2,k+q)wi,kl'l,k} 
+ {(W2,k - V2,k)('"l,k+q - ^^l,k+q) + (W2,k+q " ^^2,k+q)(wi,k " t^l.k)} (Wl,k^^l,k+q + Wl.k+q^l.k) 



85a Wo^c 



2(Ul,k?^l,k+q + Wl,k+q^^l,k)^ - Wl,k^^l,k(Wi k+q + ^l,k+q) " Wl.k+q^^l.k+ql^i k + ^^l.k) 

(Ml,k - Wl,k)(u2,k - f2,k)wi,k+qSl,k+q + (wi,k+q - Wl,k+q) (u2,k+q - W2,k+q)wi,kf l,k 



+ (W2,k - W2,k)(Ul,k+q^l,k - ^^l,k+q"l,k) + (W2,k+q - t^2,k+q) (Wl,k'*'l,k+q - ^^i^k^^l.k+q) 

, 2 

Wl,kWl,k+q + Wl,k+qf l,k ' 



(W2,k - W2,k)^Wi,k+q?^l,k+q + (W2,k+q " ?^2,k+q)^Wi,kWl,k 

+ (M2,k - ?^2,k)(u2,k+q - S2,k+q)(wi,kWl,k+q + Wl^kUl.k+q) • 



(Al) 



The corresponding vertices for the transitions within the lower (upper) band are obtained from the above expressions 
by replacing u,v ^ u,v {u,v ^ u, v) and multiplying by a factor 1/2. 
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